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Abstract 

The forward electroproduction of two light vector mesons is the first example of a colli- 
sion process between strongly interacting colorless particles for which the amplitude can 
be written completely within perturbative QCD in the Regge limit with next-to-leading 
accuracy. In a previous paper we have given a numerical determination of the amplitude 
in the case of equal photon virtualities by using a definite representation for the ampli- 
tude and a definite optimization method for the perturbative series. Here we estimate the 
systematic uncertainty of our previous determination, by considering a different repre- 
sentation of the amplitude and different optimization methods of the perturbative series. 
Moreover, we compare our result for the differential cross section at the minimum \t\ with 
a different approach, based on collinear kernel improvement. 



1 Introduction 



The BFKL approach [1] to strong interactions is expected to describe collision processes 
with a large center-of-mass energy and with a "hard" enough scale to permit the use of 
perturbative expansion in the strong coupling a s . In this approach, both in the leading log- 
arithmic approximation (LLA), which means resummation of leading energy logarithms, 
all terms (a s \n(s)) n , and in the next-to-leading approximation (NLA), which means re- 
summation of all terms a s (a s ln(s)) n , the (imaginary part of the) amplitude for a large-s 
hard collision process can be written as the convolution of the Green's function of two 
interacting Reggeized gluons with the impact factors of the colliding particles (see, for 
example, Fig. 1). 

The Green's function is determined through the BFKL equation. The NLA singlet 
kernel of the BFKL equation has been achieved in the forward case [2], after the long 
program of calculation of the NLA corrections [3] (for a review, see Ref. [4]). For the non- 
forward case the ingredients to the NLA BFKL kernel are known since a few years for the 
color octet representation in the t-channel [5] . This color representation is very important 
for the check of consistency of the s-channel unitarity with the gluon Reggeization, i.e. 
for the "bootstrap" [6]. Recently the non-forward NLA BFKL kernel has been derived 
also in the singlet color representation, i.e. in the Pomeron channel, relevant for physical 
applications [7]. 

On the side of impact factors, however, only a limited knowledge is available. Impact 
factors have been calculated with NLA accuracy for colliding partons [8] and for forward 
jet production [9]. The most important impact factor for the BFKL phenomenology, the 
7* — > 7* impact factor, is calling for a rather long calculation, although it seems to be 
close to completion now [10]. The only available colorless impact factor is presently the 
one for the forward transition from a virtual photon 7* to a light neutral vector meson 
V = p°,u,(f), obtained in Ref. [11]. This impact factor can be used, together with the 
NLA BFKL Green's function, to build completely within perturbative QCD and with 
NLA accuracy the amplitude of the 7*7* — > VV reaction. This amplitude provides us 
with an ideal theoretical laboratory for the investigation of several open questions in the 
BFKL approach and for the comparison with different approaches. 

Such investigation was started in Ref. [12], where it was first of all shown how the 
7* — > V impact factors and the BFKL Green's function can be put together to build 
up the NLA forward amplitude of the 7*7* — > VV process in the MS scheme and that 
a convenient series representation for this amplitude can be determined. Then, in the 
case of equal photon virtualities, i.e. in the so-called "pure" BFKL regime, a numerical 
study was carried out which led to conclude that the NLA corrections are large and of 
opposite sign with respect to the leading order and that they are dominated, at the lower 
energies, by the NLA correction from impact factors. However, this fact did not prevent 
us from achieving a smooth behavior of the (imaginary part of the) amplitude with the 
energy, by optimizing the choice of the energy scale sq in the BFKL approach and the 
renormalization scale which appear in subleading terms. The optimization method 
adopted there was an adaptation of the "principle of minimum sensitivity" (PMS) [13] to 
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the case where two energy parameters are present. 

The aim of this paper is to estimate the amount of the systematic effects in the deter- 
mination of Ref. [12], by exploring the two main sources: the choice of the representation 
of the amplitude and the choice of the optimization method. Concerning the first effect, in 
this paper we compare the series representation of the amplitude with another representa- 
tion, equivalent to the previous with NLA accuracy, where almost all the NLA corrections 
coming from the kernel are exponentiated. As for the second effect, we compare here the 
PMS optimization method with two other well-known methods of optimization of the 
perturbative series, namely the fast apparent convergence (FAC) method [14] and the 
Brodsky-Lepage-Mackenzie (BLM) method [15]. 

It would be quite interesting to apply to the amplitude under consideration here the 
improvement of the NLA BFKL kernel as a consequence of the analysis of collinear sin- 
gularities of the NLA corrections and by the account of further collinear terms beyond 
NLA [16, 17, 18, 19, 20, 21]. As pointed out in Ref. [12], the strategy of collinear im- 
provement has something in common with ours, in the sense that it is also inspired by 
renormalizat ion-group invariance and it also leads to the addition of terms beyond the 
NLA. Work is in progress in this direction. 

In this paper, however, we present a comparison with an approach by another research 
group based on the collinear improvement of the kernel. In particular, we compare our 
determinations for the differential cross section at the minimum \t\ of the 7*7* — > VV 
process for two values of the common photon virtuality with the results of Ref. [22] , where 
the same process has been considered using some version of a collinearly improved kernel. 

This paper is organized as follows: in the next Section we briefly recall the relevant 
notation and give the two representations of the amplitude, series and "exponentiated", 
to be considered in the following; in Section 3 we recall the strategies of the three op- 
timization methods considered and perform the numerical comparisons; in Section 4 we 
compare our differential cross section with that of Ref. [22]; in Section 5 we draw our 
conclusions. 



2 Representations of the NLA amplitude 

The process under consideration is the production of two light vector mesons (V = 
p°,u,(p) in the collision of two virtual photons, 

Y( p ) 7 *(p') - v( Pl ) v( P2 ) . (l) 

Here, neglecting the meson mass m v , pi and p 2 are taken as Sudakov vectors satisfying 
p\ = p\ = and 2{p 1 p 2 ) = s; the virtual photon momenta are instead 

p = api p 2 , p = a p 2 r pi , (2) 

as a's 
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Figure 1: Schematic representation of the amplitude for the 7*(p)7*(p') — > V(pi) V^(p 2 ) 
scattering. 

so that the photon virtualities turn to be p 2 = —Q\ and (p') 2 = —Q\. We consider the 
kinematics when 

s > Q\ 2 > A 2 QCD , (3) 

and 

a = 1 + 51 + G( s - 2 ) , a' = 1 + 51 + G( s - 2 ) . (4) 
s s 

In this case vector mesons are produced by longitudinally polarized photons in the lon- 
gitudinally polarized state [11]. Other helicity amplitudes are power suppressed, with a 
suppression factor ~ rn v /Q 12 . We will discuss here the amplitude of the forward scat- 
tering, i.e. when the transverse momenta of produced V mesons are zero or when the 
variable t — (pi — p) 2 takes its maximal value to = —Q\Q\js + (D(s~ 2 ). 

The forward amplitude in the BFKL approach may be presented as follows 

j2 -> j2-> 5+ioo 

Xm s (^) = -^/4!*i(9i,So) /^* 2 (-&ao) / (-) G u (&, q 2 ) . (5) 

5— ioo 

This representation for the amplitude is valid with NLA accuracy. Here <$>i(q\,so) and 
&2(— <?2, s ) are the impact factors describing the transitions 7*(p) — > V(pi) and 7*(p') — > 
V(p2); respectively. The Green's function in (5) obeys the BFKL equation 

<5 2 (<?i =wG f w (gi,&) - /" d 2 qK(qx,q)G u (q,q 2 ) , (6) 

where K(q\, q 2 ) is the BFKL kernel. The scale so is artificial. It is introduced in the BFKL 
approach at the time to perform the Mellin transform from the s-space to the complex 
angular momentum plane and must disappear in the full expression for the amplitude at 
each fixed order of approximation. Using the result for the meson NLA impact factor 
such cancellation was demonstrated explicitly in Ref. [11] for the process in question. 
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The impact factors are also presented as an expansion in a s 



4ire q f v 



Ms) = «.^R3(? J )+«.cSr) > D w = -w7r L v N t- 1 > W 

where fy is the meson dimensional coupling constant (f p ~ 200 MeV) and e q should be 
replaced by e/y/2, e/(3y/2) and — e/3 for the case of p°, and meson production, 
respectively. 

In the collinear factorization approach the meson transition impact factor is given as 
a convolution of the hard scattering amplitude for the production of a collinear quark- 
antiquark pair with the meson distribution amplitude (DA). The integration variable in 
this convolution is the fraction z of the meson momentum carried by the quark (z = 1 — z 
is the momentum fraction carried by the antiquark): 



The NLA correction to the hard scattering amplitude, for a photon with virtuality 
equal to Q 2 , is defined as follows 



(9) 



with t(z) given in the Eq. (75) of Ref. [11]. Ci~ 2 (q 2 ) are given by the previous expression 
with Q 2 replaced everywhere in the integrand by Q\ and Q 2 ,, respectively. We use the 
distribution amplitude in its asymptotic form ^{z) = 6z(l — z). The main reason for 
this choice is simplicity. The other point is that for the case of equal photon virtualities 
the typical values of the Reggeon momenta are q 2 ~ Q 2 . In this case the integrands in 
Eqs. (8) and (9) are smooth functions of z and, consequently, the amplitude is not very 
sensitive to the shape of the meson DA. 

In Ref. [12] the NLA forward amplitude has been written as a spectral decomposition 
on the basis of eigenfunctions of the LLA BFKL kernel: 



Xm s (A) 



(2nf 



f dv (— 

J \S 



o^{^r)ci{v)c 2 {v) 



ci{y) c 2 {v) 



+ ^)ln(^)(^) + ^) 



10 dln(^4) 



Let us briefly recall the definition of the objects entering this expression: 



at* = 



a s N c 



with N c the number of colors; 

X {y) = 2^(1) - V Q + iv) - i= Q - iv) , 



(10) 



;n) 



(12) 



Cl (i/) = Jd 2 qC[°\q 



>2\iv— 4 /-»2\ — iv— 4 



and similar equations for c^ 1 \v) and c^ 1 ^) from the NLA corrections to the impact 
factors, C ( l\q 2 ) and C^{q 2 )\ 



X{v) 



TV 2 -A 



x{v ) _ 6C (3) _ jf( v ) 



7Y 



+ 



7T 2 sinh(7rz/) 
2 1/ cosh 2 (7rz/) 



3+ 1 + ^ 



COsh(7TI/) 

71/ \ ll + 12z/ 2X 



A" c 3 / 16(1 + z/ 2 ) 



+ 4 



<j>(u) 



2 / 



cos(i/ln(a;)) 
(l + x)y/x 



7Y 



Liofx) 

6 v ; 



Li 2 (a;' 



-/ 



dt 



ln(l - f) 



(14) 
(15) 



Using Eq. (10) we construct the following representation for the amplitude 

Q1Q2 2^^s(-4.scries) 1 

I).D> s 



X 



b + ^2 ®s{hrY K hi — + d n (s , (Mr) In — 

n =l ^ \ S 0S \ s 



where the coefficients 



+00 



dv C\{y)ci{y 



Q1Q2 J ' n\ 

—00 

are determined by the kernel and the impact factors in LLA. The coefficients 

.s () \ > {) I (>„- 1 , / //]- \ n(n - 1) 



dn = nln \Q^ 



AN, 



+00 



+9l9l I du(n+l)f(u) Cl (u)c 2 (u) 



b n \Q1Q2J 2 



(n-1)! 



+00 



H ; I / dv c x {v)c 2 {v)- 



(n-l)\ 



cAv) 



c 2 {y) 



(16) 



(17) 



(18) 



are determined by the NLA corrections to the kernel and to the impact factors. Here cj^ 
are determined according to the definitions (13) from 



C^\q 2 )=C^\q 2 ) 



(19) 



1 



dz 



q 2 



q 2 + zzQ 2 



1 1 I s o \ , I (a + zz 



a 2 z 2 z 2 



— ln(a) 



Moreover, we use the notation 
5 



o { 3 \ X 3 

/(i/) = - + ^(3 + 2ii/) + ^(3 - 2iu) - V + iv ) ~ ^ (2 - ^ 



(20) 
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One should stress that both representations of the amplitude (16) and (10) are equiv- 
alent with NLA accuracy, since they differ only by next-to-NLA (NNLA) terms. The 
series representation (16) is a natural choice, since it includes in some sense the minimal 
amount of NNLA contributions; moreover, its form is the closest one to the initial goal of 
the BFKL approach, i.e. to sum selected contributions in the perturbative series. 

Actually there exist infinitely many possibilities to write a NLA amplitude. The other 
possibility considered here is to exponentiate the bulk of the kernel NLA corrections 

Tm s (A cxp ) s 7°, / s \^(^)xH+^(^)(xM+^xH[-xH+f]) 

dv — a a {ii R )ci{y)c2{y) 



x 



D±D 2 (2tt) 2 J W 

— oo 

(21) 

This form of the NLA amplitude was used in [23] (see also [24]), without account of the 
last two terms in the second line of (21), for the analysis of the total 7*7* cross section. 
We will refer in the following to this representation simply as "exponentiated" amplitude. 

It is easily seen from Eqs. (16)-(20) that the amplitude is independent in the NLA 
from the choice of energy and strong coupling scales [12]. 

For the purposes of our analysis of systematic effects, it could be acceptable also to use 
an amplitude where the last two terms in the squared bracket of the integrand in the R.H.S. 
of Eq. (21) are exponentiated. Indeed, we performed a numerical analysis also in this case 
obtaining results which are in fair agreement with our findings below. Nevertheless, we 
prefer not to exponentiate these terms for the following reason. By exponentiation we 
mean to transfer an effect from the kernel to the Green's function, which means to account 
for the corresponding effect "to all orders". The last term in Eq. (21) originates from the 
scale non-invariant part of the NLA kernel, i.e. from the running of the coupling, which 
leads in the Mellin space to the derivative term, see Eq. (25) in Ref. [12]. It is known [25], 
that an exact account of the derivative term in the BFKL equation leads to a radical 
change of both the spectrum and the eigenfunctions of the NLA BFKL kernel. Another 
point is that BFKL approach itself (where an amplitude is a convolution of the Green's 
function and the impact factors) is valid only within NLA. In higher orders one should 
take into account additional contributions, related, for instance, with the transition of 
two to four Reggeized gluons propagating in the t-channel. In this situation we decided 
to leave the last two terms in Eq. (21) unexponentiated. 

Another interesting representation for the BFKL NLA amplitude appears if one treats 
the energy scale parameters dynamically, allowing s being a function of the Reggeons 
transverse momenta. Such change of energy scale leads to the corresponding modification 
of the impact factors, and even of the kernel of BFKL equation in the case if the energy 
scale does not factorize as the product of two functions of q\ and q 2 , see [26]. How- 
ever numerical implementation of such dynamical energy scale scheme would require the 
knowledge of NLA BFKL Green's function in the momentum representation, G(q\,q2), 
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Figure 2: Im s (A ser \ es )Q 2 / {s DiD 2 ) as a function of Y for optimal choice of the energy 
parameters Yq and /j,r (curve labeled by "NLA"). The other curves represent the LLA 
result for Y = 2.2 and /ir = 10Q and the Born (2-gluon exchange) limit for /ir — Q and 
/iR = 10Q. The photon virtuality Q 2 has been fixed to 24 GeV 2 (rif — 5). 

and its subsequent integration with the impact factors. Such a study could be done with 
the methods developed in Ref. [27]. 




3 Numerical results 



In Ref. [12] we have presented some numerical results for the amplitude given in Eq. (16) 
for the Qi = Q 2 = Q kinematics, i.e. in the "pure" BFKL regime. We truncated the 
series in the R.H.S. of Eq. (16) to n — 20, after having verified that this procedure gives 
a very good approximation of the infinite sum for the Y values Y < 10 and used the 
two-loop running coupling corresponding to the value a s (Mz) = 0.12. 

We obtained there the following results for the b n and d n coefficients for rif = 5 and 
so = Q 2 = lA- 

b = 17.0664 61 = 34.5920 b 2 = 40.7609 63 = 33.0618 h = 20.7467 
b 5 = 10.5698 b 6 = 4.54792 b 7 = 1.69128 b s = 0.554475 



(22) 



d x = -3.71087 d 2 = -11.3057 d 3 = -23.3879 d A = -39.1123 
d 5 = -59.207 d 6 = -83.0365 d 7 = -111.151 d 8 = -143.06 , 



the main contribution to the d n coefficients originating from the NLA corrections to the 
impact factors for n < 3. 
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Figure 3: Xm s (A)Q 2 / (s DiD 2 ) as a function of Y at Q 2 =24 GeV 2 (rif = 5) from series 
and "exponentiated" representations, in both cases with the PMS optimization method. 

These numbers make visible the effect of the NLA corrections: the d n coefficients are 
negative and increasingly large in absolute values as the perturbative order increases. In 
such a situation it becomes necessary to optimize the perturbative expansion, by proper 
choice of the renormalization scale fi R and of the energy scale sq. 

Several ways are known to optimize a perturbative expansion. In Ref. [12] we adopted 
the principle of minimal sensitivity (PMS) [13]. Usually PMS is used to fix the value of 
the renormalization scale for the strong coupling. We used this principle in a broader 
sense, requiring the minimal sensitivity of the predictions to the change of both the 
renormalization and the energy scales, fj, R and s . In Ref. [12] we considered the amplitude 
for Q 2 =24 GeV 2 and rif — 5 and studied its sensitivity to variation of the parameters /j,r 
and Yq = ln(s /Q 2 ). We could see that for each value of Y — ln(s/Q 2 ) there are quite 
large regions in (j, R and Y where the amplitude is practically independent on ji R and Y 
and we got for the amplitude a smooth behavior in Y, shown in Fig. 2. The optimal 
values turned out to be fi R ~ 10Q and Y ~ 2, quite far from the kinematical values 
/jLr — Q and Y = 0. These "unnatural" values are a manifestation of the nature of the 
BFKL series: NLA corrections are large and then, necessarily, since the exact amplitude 
should be renorm- and energy scale invariant, the NNLA terms should be large and of 
the opposite sign with respect to the NLA. These large NNLA corrections are mimicked 
by the "unnatural" optimal values of fi R and Y . 

As anticipated in the Introduction, it is important to have an estimate of the systematic 
uncertainty which plagues our determination of the energy behavior of the amplitude. 
The main sources of systematic effects are given by the choice of the representation of 
the amplitude and by the optimization method adopted. In the following, we compare 
the determination of the amplitude at Q 2 = 24 GeV 2 (rif — 5) through the PMS method, 
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given in Fig. 2, with other determinations obtained changing either the representation of 
the amplitude or the optimization method. 

At first, we compare the series and the "exponentiated" determinations using in both 
case the PMS method. The procedure we followed to determine the energy behavior of the 
"exponentiated" amplitude is straightforward: for each fixed value of Y we determined the 
optimal choice of the parameters Hr and Y for which the amplitude given in Eq. (21) is 
the least sensitive to their variation. Also in this case we could see wide regions of stability 
of the amplitude in the (/zr, Y ) plane. The optimal values of fi R and Y are quite similar 
to those obtained in the case of the series representation, with only a slight decrease of the 
optimal hr. In Fig. 3 we show the result and compare it to the PMS determination from 
the series representation. The two curves are in good agreement at the lower energies, 
the deviation increasing for large values of Y. It should be stressed, however, that the 
applicability domain of the BFKL approach is determined by the condition a s (fj,R)Y ~ 1 
and, for Q 2 =24 GeV 2 and for the typical optimal values of fiR, one gets from this condition 
Y ~ 5. Around this value the discrepancy between the two determinations is within a 
few percent. 

We repeated the same analysis at Q 2 = 5 GeV 2 (n/ = 4) and compared the PMS 
determination from the "exponentiated" amplitude with the PMS determination from 
the series representation, obtained first in our paper Ref. [12]. Fig. 4 shows that the two 
determinations are in nice agreement. Despite that one should stress that in a case Q 2 
= 5 GeV 2 we found for exponentiated amplitude much higher value for optimal energy 
scale, Yq ~ 6. This may be an indication that the convergence of NLA BFKL approach 
is actually worse for this smaller scale than it is for Q 2 = 24 GeV 2 . 

As a second check, we changed the optimization method and applied it both to the 
series and to the "exponentiated" representation. The method considered is the fast 
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Figure 5: lm s (A scr i CS )Q 2 / (s DiD 2 ) as a function of Y at Q 2 =2A GeV 2 (n/ = 5) from the 
series representation with PMS and FAC optimization methods. 

apparent convergence (FAC) method [14], whose strategy, when applied to a usual per- 
turbative expansion, is to fix the renormalization scale to the value for which the highest 
order correction term is exactly zero. In our case, the application of the FAC method 
requires an adaptation, for two reasons: the first is that we have two energy parameters 
in the game, fiR and Y"o, the second is that, if only strict NLA corrections are taken, the 
amplitude does not depend at all on these parameters. 

Therefore, in the case of the series representation, Eq. (16), we choose to put to zero 
the sum 



and found for each fixed Y the values of fin and Y for which the vanishing occurs. 
This gives a line of values in the (/xr, Y ) plane, among which the optimal choice is done 
applying a minimum sensitivity criterion. The result is shown in Fig. 5. The agreement 
with the series representation with the PMS method is rather good over a wide energy 
range. 

In the case of the "exponentiated" amplitude", Eq. (21), we proceeded in the same 
way, but requiring the vanishing of the expression given by the R.H.S. of Eq. (21) minus 
the LLA amplitude, i.e. 



In Fig. 6 the result is compared with series representation in the PMS method: there is 
nice agreement over the whole energy range considered. 
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Figure 6: lm s (A)Q 2 /(s DiD 2 ) as a function of Y at Q 2 =2A GeV 2 (n/ = 5) from the 
series representation with the PMS optimization method and from the "exponentiated" 
representation with the FAC optimization method. 



Another popular optimization procedure is the Brodsky-Lepage-Mackenzie (BLM) 
method [15], which amounts to perform a finite renormalization to a physical scheme 
and then choose the renormalization scale in order to remove the /^-dependent part. We 
applied this method only to the series representation, Eq. (16), and proceeded as follows: 
we first performed a finite renormalization to the momentum (MOM) scheme with £ = 
(see Ref. [23]), 



a. 



l+T MOM ^ = 0) 



a. 



TT 



T, 



MOM 



(* = 0) 



rpCOnf , rp/3 

1 MOM 1 l 



MOM ? 



rpCOnf 

1 MOM 



Nc}7 1 
8 2 



T /3 

1 MOM 



2 



2.3439 



then, we chose Y and /ir in order to make the term proportional to (3 in the resulting 
amplitude vanish. We observe that the /^-dependence in the series representation of the 
amplitude is hidden into the d n coefficients, Eq. (18). Among the resulting pairs of values 
for Yq and fiR, we determined the optimal one according to minimum sensitivity. This 
method has a drawback in our case, since for each fixed Y, the optimal choice for Yq 
turned to be always Yq ~ Y . However, if one blindly applies the procedure above, one 
gets a curve which slightly overshoots the one for the series representation with the PMS 
method, see Fig. 7. 
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Figure 7: lm s (A)Q 2 /{s D X D 2 ) as a function of Y at Q 2 =24 GeV 2 (n f = 5) from the 
series representation with PMS and BLM optimization methods. 

4 The differential cross section at the minimum \t\: 
comparison with a different approach 

The 7*7* — > pp process at the lowest order (two-gluon exchange in the t-channel) was 
studied in Ref. [28]. At that level our results coincide, see also [12]. The same process 
with the inclusion of NLA BFKL effects has been considered in Ref. [22]. In that pa- 
per, the amplitude has been built with the following ingredients: leading-order impact 
factors for the 7* — > p transition, BLM scale fixing for the running of the coupling in 
the prefactor of the amplitude (the BLM scale is found using the NLA 7* — * p impact 
factor calculated in Ref. [11]) and renormalization-group-resummed BFKL kernel, with 
resummation performed on the LLA BFKL kernel at fixed coupling [29]. In Ref. [22] the 
behavior of da/dt at t — t was determined as a function of a/s for three values of the 
common photon virtuality, Q=2, 3 and 4 GeV. 

In order to make a comparison with the findings of Ref. [22], we computed da/dt at 
* = t for Q=2 and Q=4 GeV as functions of y/s. We used / p =216 MeV, a EM = 1/137 
and the two-loop running strong coupling corresponding to the value a s (M z ) = 0.12. The 
results are shown in the linear-log plots of Figs. 8 and 9, which show a large disagreement. 
It would be interesting to understand to what extent this disagreement is due to the use 
in Ref. [22] of LLA impact factors instead of the NLA ones or to the way the collinear 
improvement of the kernel is performed. 

In order to understand to what extent the discrepancy is due to the use of leading order 
(LO) impact factors instead of next-to-leading order (NLO) ones, we repeated our deter- 
mination of da/dt at t = t for Q=2 and Q=A GeV, using LO impact factors and keeping 
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Figure 8: Linear-log plot of da / dt\ t=tQ [pb/GeV 2 ] as a function of ^/s at Q 2 =1Q GeV 2 
{rif — 4) from the series representation with the PMS optimization method (solid line) 
compared with the determination from the approach in Ref. [22] (dashed line). 

from the their NLO contribution only the terms proportional to ln[s /(<5i<52)] and to 
\n[fi 2 R /(QiQ 2 )] which are universal and needed to guarantee the s - and /^-independence 
of the amplitude with NLA accuracy. The result is that da/dt at t = t Q increases roughly 
by an order of magnitude with respect to our previous determination (see Figs. 10 and 11) 
and therefore the disagreement with [22] becomes even worse. This is not surprising: im- 
pact factors give a sizable contribution to the NLA part of the amplitude which is negative 
with respect to the LLA part; if they are kept at LO, the NLA part of the amplitude is 
less negative and the total amplitude is therefore increased. 



5 Conclusions 



We have studied the amplitude for the forward transition from two virtual photons with 
equal virtuality to two light vector mesons in the Regge limit of QCD with next-to-leading 
order accuracy. We have found that its behavior with the center-of-mass energy is stable 
in the applicability region of the BFKL approach under change of representation of the 
amplitude and under change of the method of optimization of the perturbative series. 

We have determined also the differential cross section at the minimum \t\ for two 
values of the common photon virtuality and found strong disagreement with another 
determination based on the inclusion of next-to-leading order effects only through the 
kernel and on the use of collinear improvement of the kernel at the leading order. 
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Figure 9: Linear-log plot of da/dt\ t=t [pb/GeV 2 ] as a function of ^/s at Q 2 =A GeV 2 
(rif — 3) from the series representation with the PMS optimization method (solid line) 
compared with the determination from the approach in Ref. [22] (dashed line). 
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